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Abstract 

By using the quantum kinetic approach with the instantaneous local equilibrium approximation, 
we propose an equation that is capable of addressing magnetization dynamics for a wide range of 
temperatures. The equation reduces to the Landau-Lifshitz equation at low temperatures and to 
the paramagnetic Bloch equation at high temperatures. Near the Curie temperature, the magneti- 
zation reversal and dynamics depend on both transverse and longitudinal relaxations. We further 
include the stochastic fields in the dynamic equation in order to take into account fluctuation at 
high temperatures. Our proposed equation may be broadly used for modeling laser pump-probe 
experiments and heat assisted magnetic recording. 

PACS numbers: 75.78.-n, 75.40. Gb 



I. INTRODUCTION 



The phenomenological Landau-Lifshitz (LL) equation is the basis of powerful micromag- 
netic codes for simulation of magnetic structure and dynamics in magnetic materials. The 
key component of the LL equation is that magnetization relaxation during dynamic processes 
is described by a single damping parameter a PQ, 



where m(r, t) is the magnetization density vector which is a function of space and time, 
m = |m| is its magnitude, and H e g is the effective magnetic field including the magnetic 
anisotropic, magnetostatic and external fields. The second term on the right side of the 
equation describes a phenomenological transverse relaxation since the magnitude of the 
magnetization density m is conserved. Such transverse relaxation model is indeed a valid 
approximation because the magnetization m (the order parameter) of the ferromagnet is 
nearly independent of the magnetic field as long as the temperature is not too close to the 
Curie temperature. 

Recently, there is an emerging technological need to extend the LL equation to high 
temperatures in order to model the dynamics near or above the Curie temperature for laser- 
induced demagnetization (LID) [21 E] and heat-assisted magnetic recording (HAMR) (H [5] . 
Due to the strong fluctuation of the magnetic momentum, one needs to reduce the size of 
magnetic cells if one continues to use the LL equation to model the magnetization dynamics. 
When the cell size reduces to the ultimate smallest size of magnetic atoms, the so-called 
atomistic LL equation which has the same form as the conventional LL equation [6] had 
been proposed, 



where H is an effective magnetic field (treated as a c-number) including a random fluctuating 
field and Sj is the spin of ith atom which is treated classically. While the above atomistic LL 
equation might qualitatively capture some of static and dynamic properties near the Curie 
temperature [7J E], we point out below that the above atomistic LL equation has several 
fundamental problems. 

First, the spin Si in ferromagnetic metals such as Ni, Co, Fe and their alloys is usually 
small. The replacement of the quantum spin by the classical vector severely neglects the 
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quantum nature of the spin fluctuation of atomic spins. More importantly, the atomistic 
LL equation, Eq. (2), has no microscopic origin and it is fundamentally incompatible with 
quantum mechanics. For example, if one takes the case for Si — 1/2 (e.g., Ni). The second 
("damping") term of Eq. (2) becomes 7«(|SjXH — |H) for a quantum spin and thus Eq. (2) 
becomes 

-^ = -7(l + ia/2)SiXH + 7aH/2. (3) 

This unphysical equation originates from the broken time-reversal symmetry inherited on 
the atomistic LL equation. 

The second difficulty is that the atomistic LL equation is not derivable from an effective 
Hamiltonian, even at the phenomenological level. If the atomistic LL equation has some 
validity, a microscopic or an effective Hamiltonian should exist. For example, if we construct 
a spin Hamiltonian of the form H' oc ^ Sj • (H + aSj x H), the equation of motion for Sj 
would be dSi/dt = (l/ih)[Si, H'] which results in an additional term compared to Eq. (2) 
due to none-zero commutation [Sj, aSj x H] ^ 0. On the other hand, if one takes the 
phenomenological Hamiltonian as H' cx £\ Si • (H + am x H) where m =< Sj > (note that 
<> denotes ensemble thermal averaging and thus m is a c- number), the result dynamics for 
Sj would be 

^ = i[Si, H'} = - 7 S, x H - 7 «Sj x (m x H). (4) 

The above equation is precisely the original LL equation after the thermal averaging. Thus, 
the macroscopic LL equation is derivable from an effective Hamiltonian while the atomistic 
LL equation is not. 

In spite of above conceptual difficulties in the atomistic LL equation, it has been shown 
that the result derived from the atomistic LL equation with stochastic fields is in agreement 
with the Monte Carlo simulation [9]. We point out that this agreement is not surprising: 
for equilibrium properties such as magnetic moment and critical exponents, the calculated 
results are insensitive to the details of the "damping" ; for dynamic properties such as reversal 
time, the Monte Carlo steps are calibrated to fit the real time in the atomistic stochastic 
LL equation (TO]- Therefore, such agreement should not be interpreted as the proof of the 
validity of the atomistic LL equation. 

In this paper, we propose an effective magnetization dynamic equation for a wide range 
of temperatures without assuming the presence of the atomistic LL equation for each atomic 
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spin. By using the equation of motion for the quantum density matrix within the instanta- 
neous local relaxation time approximation [TT], we show that the magnetization dynamics 
for ferromagnets can be cast in the form of the Bloch equation for paramagnetic spins [T2] . 
In Sec. II, we explicitly derive the generalized Bloch equation and show that the equation is 
consistent with the known dynamics at low and high temperatures. In Sec. Ill, we analyze 
the longitudinal and transverse relaxations from our result, and apply our effective equation 
to study the magnetization reversal processes near Curie temperatures. Finally, we add 
necessary stochastic fields in the equation to capture the fluctuation of the dynamics. 



II. EFFECTIVE DYNAMIC EQUATION FOR FERROMAGNETS 

We start with a density operator p which may be written in the spinor form p = pi + cr- p 2 
where p\ and p 2 are spin-independent and spin-dependent density operators, and cr is the 
Pauli matrix vector. Within the instantaneous local relaxation time approximation, the 
density operator satisfies the quantum kinetic equation [TT] 

d P 1 r- Ai Pi - Pi P2~ P2 / r \ 

~dt = ^ih Tp ' r\j (5) 

where p\ and p 2 are the instantaneous local equilibrium (ILE) densities; they are different 

from the static equilibrium values. In electron transport theories, these ILE densities depend 

on the local chemical potential p(r) or the local electric field E(r,t) and they are in turn 

related to the densities themselves. For example, for spin dependent electron transport, the 

inclusion of the spin relaxation (third term of Eq. 5) leads to the well-known spin-diffusion 

equation for the spin dependent chemical potential (or spin density) [13]. In the present 

case, these ILE densities are functions of the local effective magnetic field. At a given time, 

the effective field consists of the ferromagnetic exchange, anisotropy, external, and classical 

magnetostatic field; we will discuss these fields in more details later. The two relaxation 

times t p and t s represent the momentum and spin relaxation times; these two relaxation 

times control the electron charge diffusion (conductance) and spin diffusion (spin-dependent 

transport). If we now consider an effective Hamiltonian H = H$ — gp,cr ■ H t (t) where p is the 

Bohr magneton, H is treated as an unperturbed Hamiltonian, we find the self-consistent 

equation for the magnetization m = gpTi{crp) = gpTrp 2 readily from Eq. (5), 

dm m — m ra (H t ) , , 

— = - 7 mxH t 6 

at t„ 



where the ILE magnetization m eq = g^p 2 is identified as the thermal equilibrium value for 
a given magnetic field H t . 

At the first sight, Eq. (6) is similar to the well known Bloch equation [T2] that has been 
widely used for understanding nuclear spin resonance experiments. In the Bloch equation, 
the equilibrium magnetization m eq is a known equilibrium state which is related to the 
dynamic susceptibility i- e -> m e<? = xH ex t and m eq is independent of m(£). In the 

present content, m eq is not known a priori and m eq varies with time. At any time t, there 
is an instantaneous equilibrium magnetization m eq that depends on the total magnetic field 
H t . To solve Eq. (6), one first needs to model the instantaneous local field H t and its relation 
to m cq . 

In the conventional LL equation, the effective field H e e consists of the external field, the 
anisotropy and the magnetostatic (dipole) fields. The exchange field which comes from the 
ferromagnetic exchange interaction between neighboring spins is included only when there 
is spatial variation in the magnetic domain structure. The uniform exchange term, Jm, is 
unimportant since it is parallel to the magnetization and it does not contribute to the LL 
dynamic equation. In the present case, however, the exchange interaction is the largest and 
most important term in determining the instantaneous equilibrium magnetization m eq . We 
thus model the total instantaneous magnetic field H t = Jm + H e g. It is noted that H c g 
depends on the instantaneous magnetization m(£) as well. 

Next, we should establish an explicit relation between the total field H t with m eq . There 
are a number of approaches available to describe such relation. The simplest approach 
would be using the molecular field approximation where the equilibrium magnetization can 
be explicitly expressed by [14J 

m eq = gji < Si >= gnSBs{f3gnH t )H t (7) 

where S is the spin of the atom, (3 = (fc^T) -1 is the inverse of temperature, Bs(x) = 
(1/S)[(S + l/2)coth(S + i/2)x - (l/2)coth(z/2)] is the Brillouin function and H t is the 
unit vector in the direction of H t , i.e., H t = Ht/H t . In the time-independent case, m = m eq 
and the above equation is the well-known mean-field result that determines the ferromagnetic 
order parameter m eq . In a non-equilibrium situation where m depends on time, we interpret 
m eq in Eq. (7), which is also dependent on time, as the instantaneous local equilibrium 
magnetization at a given (instantaneous) field H t . 
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Our proposed Eq. (6) supplemented by Eq. (7) can semi-quantitatively describe magneti- 
zation dynamics at all temperatures. Before we examine some limiting cases, we comment on 
certain important approximations leading to these equations. The instantaneous relaxation 
time approximation, Eq. (5), has been routinely applied to many quantum or semi-classical 
systems for transport and magnetic properties. The accuracy of this approximation is hard 
to assess for the ferromagnetic systems. However, the instantaneous relaxation time approx- 
imation has been very successfully applied in spin diffusion of magnetic multilayers where 
the semiclassical distribution function is assumed to relax to the instantaneous chemical 
potential [32]. Furthermore, the relaxation time approximation usually serves as a first step 
in a phenomenological theory since it gives rise an analytically closed form. The most severe 
approximation is to replace m cq by the mean field Brillouin function, Eq. (7). Such ap- 
proximations are known to produce inaccurate critical exponents and Curie temperatures. 
There are several much improved approaches such as the renormalization group theory [15J, 
self-consistent random phase approximation [TB], and Monte Carlo simulation [T7j. While 
these approaches treat the fluctuation near the critical temperature better, they are far more 
complicated and without an analytical form. On the other hand, the mean field approx- 
imation is qualitatively correct and it allows a much simpler description of magnetization 
dynamics in spite of underestimating the critical fluctuation. For the purpose of establishing 
a phenomenological dynamic equation similar to the LL equation, we believe that the choice 
of the mean field approximation throughout this study is appropriate. 

Similar to the LL equation, Eq. (6) contains a phenomenological parameter, r s , repre- 
senting the magnetic relaxation of paramagnetic spins. In transition metals, t s is related to 
the spin-flip time. In fact, there are a number of theoretical and experimental studies on the 
numerical values of r s in different materials [TBTI20] . For transition metals, the relaxation 
time ranges from sub-picoseconds to a few picoseconds. 

III. LONGITUDINAL AND TRANSVERSE MAGNETIZATION DYNAMICS 

Before we proceed to solve Eq. (6) in a number of interesting examples, we examine 
several limiting cases. First, by using the identity 
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we write Eq. (6) in terms of three mutually perpendicular vectors, 

dm 7a tr . , 70;; , TT . 

— = -7mxH cff m x (m x H cff ) (m • H t )m (9) 

at m m 

where we have introduced the transverse and longitudinal dimensionless damping coefficients 
CCtr an d ai, 

= (10) 

and 

f 

cti = — 

7T. 



m m eq 



m • H t mH t 

At low temperatures, m is close to gfiS and the exchange field Jm is much larger than 
the other fields H e fj. Thus, one immediately has a tr = (7T s Jm) _1 . In a typical transition 
ferromagnet such as Co or Fe, J is of the order of the Curie temperature (0.1-0.2 eV) and 
t s is a sub-picosecond, we find a tI is of the order of 10" 3 — 10 _1 . 

To estimate the low temperature longitudinal relaxation ai from Eq. (11), we consider 
an initial m deviates from the equilibrium value of gfiSB s and from Eq. (11), ol\ would be 
about the same order of magnitude as a tT . However, the longitudinal field Jm is much larger 
than H c ff and thus the ratio of the longitudinal (tj) to the transverse (r tr ) relaxation times is 
about r//r tr ~ H c s/J. Even for a very high anisotropy material and a large magnetic field, 
J is several orders of magnitude larger than H e s; this justifies that at the low temperature 
one can neglect the longitudinal relaxation in the dynamic equation, i.e., the magnitude of 
the magnetization is always in equilibrium. 

When the temperature is much higher than the Curie temperature, Eq. (6) represents 
the paramagnetic Bloch equation. In this case, the equilibrium magnetization m eq may 
be expressed via susceptibility Xi m e<? — X^-cS- Such dynamic equations have been 
frequently used for understanding paramagnetic resonant phenomena where the resonance 
width is determined by the relaxation time t s . 

The most interesting case of Eq. (6) is for temperature close to Curie temperature where 
transverse and longitudinal relaxation times could become comparable. To see this, we 
consider the effective field is parallel to m(t) = m{t)e z and expand Bg{x) — (S + l)x/3 — 
(1/90) (5 + 1)(2S 2 + 2S+ l)x 3 up to the third order in the small x where x = /3gfiH t . Then, 
Eq. (6) for temperature close to the Curie temperature becomes 



dm 



dt Jr. 



TA 3 T c 3 / 1 

T ) * 10J 2 T 3 VS 2 ' (1 + S) 



" * H t + — ^ ( — + ] Hf - H cS 



(12) 
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FIG. 1: (Color online) The longitudinal relaxation time 17 as a function of temperature for several 
magnetic fields. We choose a small difference between m and m eq at t = and identify t = r\ 
where the difference is reduced by the half. We have used t s = 1 ps and S = 1/2. 

where T c = S(S + 1) J(gfi) 2 /3k B is the mean field Curie temperature. In the absence of the 
magnetic field H eff = and H t = Jm, and we can immediately solve the above equation, 

-1/2 



m(t) = m(0)e- t/Tl [l + G(l - e~ 2t/n 
where m(0) is the initial magnetization, G 



(13) 



iot 3 



3T e 3 m 2 (0) ( J_ l_ 

S 2 "r" (1+5) 2 



1 - T c /T)-\ and 



T» 1 



7; 



(14) 



Thus, the longitudinal relaxation time, |tj|, near Curie temperature, is associated with the 
critical phenomenon. The relaxation time becomes very long when the temperature ap- 
proaches the Curie temperature. The dynamics slow-down at the critical temperature is in 
fact a general property of critical phenomena [21J. In the presence of the magnetic field, the 
phase transition becomes a smooth change and the dynamic slow-down is no more critical. 
In Fig. 1, we show the longitudinal relaxation as the function of the magnetic field and 
temperature. Clearly, the magnetic field suppresses the longitudinal dynamic slowdown. It 
is noted that the peak of the relaxation time in the presence of the magnetic field is shifted 
to higher temperatures. 

In order to gain more quantitative insight for the interplay between the transverse and 
longitudinal relaxations, we consider several simple cases where the numerical calculations 
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FIG. 2: (Color online) a) The time dependence of magnetization reversal when a reversal magnetic 
field is applied at t > for the temperature T = 0.9T C and T = 1.0T C as indicated. The blue 
(solid) and red (dashed) curves were obtained by Eq. (6) and by the first two terms of Eq. (9), 
respectively, b) the reversal times and their ratio as a function of the temperature obtained from 
a). The parameters are t s = 1.0 ps, S = 1/2, the anisotropy constant K = 0, and the external 
field H = 1.0 T. 

can be readily performed. We assume that the magnetic particle is a single domain so that 
there is no spatial dependence of the effective field and the magnetization. Furthermore, the 
long-range magnetostatic field is also discarded. In the first case, we compare the reversal 
times with and without the longitudinal relaxation in a simplest case: an isotropic magnetic 
particle (zero magnetic anisotropy) is initially magnetized at 5° from +z axis and a reversal 
magnetic field in the direction of —z is applied at t > 0. Figure 2(a) shows the importance 
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of the longitudinal relaxation when the temperature approaches Curie temperature. We 
compare the magnetization dynamics with and without the last term of Eq. (9). If the 
temperature is considerably below the Curie temperature, e.g., T = 0.9T C , the longitudinal 
relaxation term has a negligible effect, i.e., the result is essentially same whether the last term 
of Eq. (9) is included. This is because the magnitude of the magnetization is nearly time- 
independent at low temperature. When the temperature is near the Curie temperature, the 
magnitude of the magnetization is significantly reduced. More importantly, the magnitude 
is now a function of time due to its dependence on the total effective field. In this case, there 
is a much difference if one includes the longitudinal relaxation. In Fig. 2(b), we show the 
ratio of the reversal times calculated with and without the longitudinal relaxation. Clearly, 
the reversal time from Eq. (6) is much faster than that of the LL equation if the temperature 
is close to or higher than Curie temperature. 

Next we apply our equation to a hypothetical HAMR process when the laser heating and 
thermal diffusion produce a time-dependent temperature profile: the temperature of the 
particle increases linearly T(t) = T rm + (t/t heat )(T p — T rm ) from the room temperature T rm to 
a peak value T p for the period of < t < theat of lasing application. After the heating process 
is completed and the laser is removed, the temperature decreases due to heat diffusion into 
surroundings. We assume the temperature is T(t) = T rm + (T p — T rm ) exp[— (t — theat)Acooi] 
for t > thcat- While the precise temperature profile should be determined via heat transport 
equations with proper boundary conditions, our hypothetical temperature is characterized 
by three parameters: the peak temperature of the particle T p , and the heating and cooling 
rates 1/iheat and l/t CO oi- We choose the low-temperature magnetic anisotropy field much 
larger than the external magnetic field so that the magnetic reversal does not occur at the 
room temperatures. The temperature dependence of the anisotropy energy E a is modeled by 
E a = Km 2 (T) sin 2 9, where m(T) is the magnitude of the magnetization at temperature T 
and 9 is the angle between the magnetization vector and z-axis [221 123]- By placing the above 
temperature profile and effective magnetic fields into Eq. (6), we have numerically calculated 
the time dependent magnetization shown in Fig. (3). As we expected, the magnetization 
reversal requires a high peak temperature T p to reduce the anisotropy. The rates of heating 
and cooling are also important; they should be slow enough so that the magnetization has 
sufficient time to relax to the ground state via transverse and longitudinal relaxations. 

More quantitatively, we have made two comparisons in Fig. 3. First, we compare our 
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FIG. 3: (Color online) Time dependence of magnetization for given temperature profiles after a 
reversal magnetic field H = 1.5T is applied for t > 0. a) The results are obtained from Eq. (6) 
(solid curves) and from Eq. (15) (dashed curves), b) Comparison of the results from Eq. (6) and 
from the conventional LL with a constant damping parameter (dashed curves). The Inserts are 
the hypothetical temperature profiles. The parameters are K = 8.0T, S = 1/2, and t s = 1.0 ps. 



equation with a modified LL which allows the magnitude of the magnetization varying with 
the time due to changing temperature in HAMR, 



dm d(m cq m 



dm Pn dT 



"eq 



-7m x H cff - 7 — m x (m x H cff ) + 

m dT dt 



m 



(15) 



dt dt 

where the transverse damping parameter is given by Eq. (10). The above equation implies 
that the magnitude of the magnetization is always in equilibrium with the instantaneous 
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temperature, i.e., the longitudinal relaxation is infinite fast. Fig. 3(a) shows that such ap- 
proximation is quite accurate even for the temperature T p = 0.95T C . However, the deviation 
begins to show up when the peak temperature is higher than the Curie temperature. In 
Fig. 3(b), we further compare our results with a constant damping parameter (i.e., taking 
atr in Eq. (15) as a constant). The deviation of this conventional LL with ours becomes 
more significant at high temperatures. For example, even for T p = 1.2T C , the magnetization 
reversal is not possible from the conventional LL equation, see Fig. 3(b). 

To end this section, we should briefly compare our equation with the LLB equation of 
Garanin [B] . Since the LLB equation is based on the atomistic LL equation that we believe 
is questionable, we should not make extensive comparisons. We point out that the LLB 
equation also contains the transverse and longitudinal relaxations, and the essential differ- 
ence is the temperature dependence of the relaxation parameters. At low temperatures, 
both our equation and the LLB equation reduce to the conventional LL equation. At high 
temperatures, the relaxations in the LLB equation depend explicitly on temperatures; this 
is because the longitudinal relaxation to the equilibrium magnetization is solely controlled 
by the classical random field which is proportional to the temperature. In our case, the de- 
pendence of the relaxation time on temperature is implicit, via the temperature dependence 
of the equilibrium magnetization. More importantly, the instantaneous relaxation time t s 
in our equation has microscopical meaning as the scattering lifetime of electron spins while 
the damping parameter in the atomistic LL does not have a microscopic counterpart. 

IV. STOCHASTIC FIELDS 

Our proposed equation, Eq. (6), describes the time-dependence of the average magnetiza- 
tion. The fluctuation at the finite temperature, particularly at a high temperature, becomes 
important. To address the fluctuation, one should include stochastic fields in the macro- 
scopic dynamic equation. Similar to Brown's method [24J for the LL equation, we introduce 
the stochastic fields h(t) as follows, 



We point out that the stochastic field does not enter in the relaxation term although the 
instantaneous equilibrium magnetization depends on the total field H t . The reason is as 




dm 



m — m, 



(16) 
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follows. The interaction between the random field and the magnetization is — m • h(t). This 
interaction gives arise a random torque on the magnetization —7m x h(t) that is added to 
the deterministic torque equation. As in the case of the Brownian motion, the Langevin 
random field f(t) is only included in the particle motion d 2 r/dt 2 = — av + F(t) + f(t) (where 
the friction force -av and the external driven force F(i) are not changed by the random 
force). To determine the magnitude and the correlation of the stochastic fields, we first write 
the above stochastic equation in the standard form of Langevin, 

drrii 



dt 



-7m x H cff 



m — m 



cq 



JZ-jZjjk'mjhk- 

jk 



(17) 



where is the Levi-Civita symbol. The corresponding Fokker-Planck equation is thus 
have the following form, 



dP 
~dt 



d 
drrij 



—7m x H, 



m — m 



cff 



eq 



+ £>7 2 m x ( m x 



d 
dm 



(18) 



where P is the probability density and D is the random field correlation constant. At 
the equilibrium, one may assume that the probability density takes a simple Boltzmann 
distribution, i.e., P oc exp(— m • H). By placing this form of P into Eq. (18), one finds the 
desired correlation of the random field given below, 



1 um /n\ 2k B Ta tr 



(19) 



V. SUMMARY 



In this paper, we have proposed a model of magnetization dynamics for an entire range of 
temperature based on the quantum kinetic approach with the instantaneous local relaxation 
time approximation. The resulting equation generates a low temperature magnetization 
dynamic same as the Landau-Lifshitz equation, namely, the transverse magnetization is 
sufficient to describe dynamics. When the temperature approaches or exceeds the Curie 
temperature, it is essential to include the longitudinal magnetization relaxation. With our 
new dynamic equation, one can model the entire heat-assisted magnetic recording processes 
when the temperature are heated and cooled through the Curie temperature [25H27] . The 
stochastic fields on the magnetization are also proposed. This work is partially supported 



13 



by the U.S. DOE (DE-FG02-06ER46307) and by the NSF (ECCS-1127751). 



[1] L. D. Landau and E. M. Lifshitz, Phys. Z. Sowietunion 8, 153(1935). 

[2] E. Beaurepaire, J. C. Merle, A. Daunois and J. Y. Bigot, Phys. Rev. Lett. 76, 4250(1996). 
[3] Stanciu, C. D. and Hansteen, F. and Kimel, A. V. and Kirilyuk, A. and Tsukamoto, A. and 

Itoh, A. and Rasing, Phys. Rev. Lett. 99, 047601(2007). 
[4] M. H. Kryder et al., Proc. IEEE 96, 1810 (2008). 
[5] W. A. Challener, et al., Nature Photonics 3, 220 (2009). 
[6] D. A. Garanin, Phys. Rev. B 55, 3050 (1997). 

[7] O. Chubykalo-Fesenko, U. Nowak, R. W. Chantrell, and D. A. Garanin, Phys. Rev. 74, 094436 
(2006). 

[8] U. Atxitia, et al, Phys. Rev. B 82, 134440 (2010). 

[9] U. Nowak, R. W. Chantrell, and E. C. Kennedy, Phys. Rev. Lett. 84, 163 (2000). 
[10] X. Z. Cheng et al, PRL 96, 067208 (2006). 
[11] R. W. Davies and F. A. Blum, Phys. Rev. B 3, 3321 (1971) 
[12] F. Bloch, Phys. Rev. 70, 460 (1946). 
[13] T. Valet and A. Fert, Phys. Rev. B 48, 7113 (1993). 

[14] N. W. Ashroft and N. D. Mermin, Solid State Physics (Holt, Rinehart and Winston, New 
York) 1976. 

[15] K. G. Wilson, Phys. Rev. B 4, 3174 (1971); Phys. Rev. B 4, 3184 (1971). 
[16] Norberto Majlis, The Quantum Theory of Magnetism (World Scientific, Singapore) 2000. 
[17] K. Binder and D. W. Heermann, Monte Carlo Simulation in Statistical Physics (Springer- 
Verlag, Berlin) 1997. 

[18] D. V. Fedorov, P. Zahn, M. Gradhand and I. Mertig, Phys. Rev. B 77, 092406 (2008). 
[19] P. Monod and S. Schultz, J. Phys. (Paris) 43, 393 (1982). 

[20] A. C. Gossard, A. J. Heeger, and J. H. Wernick, J. Appl. Phys. 38, 1251 (1967). 
[21] P. C. Hohenberg und B. I. Halperin, Rev. Mod. Phys. 49, 435 (1977). 

[22] J. U. Thiele, K. R. Coffey, M. F. Toney, J. A. Hedstrom, and A. J. Kellock, J. Appl. Phys. 

14 



91, 6595 (2002). 

[23] P. Asselin, R. F. L. Evans, J. Barker, R. W. Chantrell, R. Yanes, O. Chubykalo-Fesenko, D. 

Hinzke, and U. Nowak, Phys. Rev. B 82, 054415 (2010). 
[24] W. F. Brown, Phys. Rev. 130, 1677 (1963). 

[25] A. Torabi, J. van Ek, E. Champion, and J. Wang, IEEE Trans. Magn. 45, 3848 (2009). 
[26] C. Bunce et al, Phys. Rev. B 81, 174428 (2010). 

[27] J. I. Mercer, M. L. Plumer, J. P. Whitehead, and J. van Ek, Appl. Phys. Lett. 98, 192508 
(2011). 



15 



